The electric double layer structure around charged spherical interfaces 
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We derive a formally simple approximate analytical solution to the Poisson-Boltzmann equation 
for the spherical system via a geometric mapping. Its regime of applicability in the parameter space 
of the spherical radius and the surface potential is determined, and its superiority over the linearized 



o 

^Nj ■ solution is demonstrated. 



Charged objects in electrolyte solutions are surrounded by electric double layers (EDL}i. One ionic layer is due to 
a host of chemical interactions, and the second layer is formed by the excess of oppositely charged ions in the solution. 



f"]^ ■ screening the charged objects. The EDL structure is responsible for the stability of colloidal dispersions^, and various 

o" 

' electrostatic phenomena in biophysical systems^. The distribution of the screening potential in EDL is characterized 



by the Poisson-Boltzmann (PB) equation. In this letter, we will study the EDL structure around charged spherical 



g ■ Dy tne roisson-eoitzmann [t^u ) equation, m tnis letter, we win st 

■ interfaces, which are ubiquitous in colloidal and biophysical systems^ 
O ■ Largely due to its nonlinear nature, analytical solutions to the PB equation are available only for planar— and 

cylindrical^ systems. For a spherical system, the analytical solution to the linearized PB equation is available^. 
^ ] Despite a recently proposed analytical series solutionS., a formally simple approximate solution is still in demand 
[ for studying analytical problems that are based on the screening potential in EDL. Numerical techniques^ and the 

■ 

Debye-Huckel linearized approximations have long been the only available basic methods to solve the PB equation. 
Various perturbative solutions have been proposed based on the planar solutioiii^, or the linearized solutionii, as the 
zeroth order approximation. Perturbative methods that start with the planar solution are limited to the regime of 

■ large spherical radius, while those which start with the linearized solution work in the weak potential regime. The 
. . ■ geometric construction of a formally simple approximate analytical solution that can match both the planar and the 

linearized solutions is one concern of this letter. 

In q:-q symmetric electrolytes, the dimensionless Poisson-Boltzmann equation is2 

Alp = s'mhij}, (1) 

in the natural units ksT/q and the Debye length = eksT/ (Simq"^). q is the absolute value of the charge of ions. 
We propose an analytic scheme to yield an approximate solution to the Poisson-Boltzmann equation in a spherical 
system from the known analytical solution to the planar problem. The planar and spherical systems are connected by 



d"("z)T 



z 



P 



z 



FIG. 1: The mapping from the z > bulk space to the space outside a sphere with radius a (represented by the blue circle) 
via. consecutive stereographic projections from a plane to a sphere. A deviation d{z) is introduced to guarantee the equality 
of the potential on the blue spherical shell and the potential at the plane located at z. 

a geometric mapping defined in Fig. [TJ Both the planar and spherical systems are schematically plotted in Fig. [T] The 
planar system is composed of a charged plate locating on the x-y plane with the solution in the z > bulk space; and 
the corresponding spherical system consists of a charged sphere of radius a centered at the origin, which is immersed 
in the solution in the r > a bulk space. The surface potential in both systems is denoted as ■00 • The potential in the 
planar system is known as ip{z), while that around the corresponding spherical system is tpir). These two potentials 
can be related by a geometric mapping which is defined as follows. By moving an arbitrary equipotential plane at z in 
the planar system by d{z) followed by a stereographic projection as shown in Fig. [l] the equipotential plane ip{z) in 
the plate system is geometrically mapped to the equipotential spherical shell Tjj{r(z)) in the corresponding spherical 
system. The charged plate on the x-y plane is mapped to the spherical interface at r = a, so d{z 0) = 0. The whole 
z > bulk space in the planar system can be mapped to the bulk space outside the spherical interface by repeating 
the mapping defined above for all equipotential planes below the x-y plane. 
The displacement field d{z) is introduced to guarantee that 



with r = a + z — d{z). d{z) encodes all information of the potential about the spherical interface. Geometrically d{z) 
describes how the equipotential planes in the planar system squeeze to form the equipotential spherical shells in the 
corresponding spherical system. The problem for solving the PB equation in a spherical system is now converted to 
solving for the geometric deviation d{z). The form of d{z) can be found in the weak potential regime where both 
4'l{z) of a planar system and of a spherical system are known: 



tp{z) = ip{r) 



(2) 



iiL{z) = V'oe 



(3) 
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FIG. 2: (a) The plot of the displacement field d{z) for various spherical radii: a — 0.1 (green), a = 1 (red), a = 10 (blue), 
a = 100 (black). As a ^ oo, d{z) is expected to vanish, (b) The plot of the strain fields d'{z) for the corresponding displacement 
fields d(z) in (a). 



and 



(4) 



where a is the radius of the spherical interface and the subscript L stands for linearized solution. Inserting Eq. ([3l |4]) 
into Eq. ([2]) leads to a/{a + z — d{z)) — cxp(— d(z)), from which we have 



d{z) = a + z-W{ae°'+''), 



(5) 



where W{x) is the Lambert's W function defined by x = W{x) exp{W{x))^. It is checked that d{z ^ 0) = and 
d{z) as a — cxD. And d{z oo) ^ ln(z/a), since asymptotically W{x — > oo) ~ \nx — ln(hia;) (Ref. [12]). 
The Lambert's W function is also found in other physical systems, such as the fringe field of a capacitor and Wiens 
displacement law in black body radiationi^. Eq. ([S]) shows that d(z) is independent of Tpo in the weak potential limit, 
since ipo appears as a prefactor in both V'(^) and in the weak potential limit as shown in Eq. ([31131). 

The plot of d(z) for various spherical radii is given in Fig. [21 (a). The squeezing of equipotential surfaces near 
a spherical interface with smaller radius is seen to be larger. It is expected that the displacement vanishes for an 
infinitely large spherical interface that approaches a plate. Fig. [51 (a) also shows the behaviors of d{z) in two regions 
- steep slope for small z and much smaller slope for large z where d{z — > oo) ~ ln(z/a). It gives the qualitative 
picture of forming equipotential spherical shells from the corresponding equipotential planar slices. Near the charged 
spherical interface each slice moves more than the slice ahead of it that is closer to the interface. In the region far 
away from the interface, the squeezing is much more uniform. The relative squeezing of equipotential surfaces is 
better demonstrated in terms of the "strain" field d'{z). The strain field plotted in Fig. [21(b) shows that the strain 
is concentrated near the spherical interface and the strain concentration is more significant near spherical interfaces 
of smaller radii. 
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FIG. 3: The potential around a charged spherical interface of unit radius in a q:q solution. The surface potentials in (a)-(d) are 
respectively 0.1, 1,5, 10 measured by q/iksT) = 1. The curves give the numerical solution to the Poisson-Boltzmann equation 
(black), analytical solution for weak ipo (green dashed) and the solution constructed by the our method (red dots). 

We suggest that the form of d{z) for arbitrary be approximated by Eq. ([5]) under the assumption of weak 
dependence of d{z) on -00 • This assumption is to be substantiated later. We can then construct the analytical 
approximate solution to the PB equation for a spherical system from the known analytical solution to a planar 
system, which is 

1 +76""^ 



■0 = 2 In ■ 



(6) 



1—76 ^ 

where 7 = (exp(?/'o/2) — l)/(exp(-0o/2) + 1). On the other hand, r[z) ~ M^(aexp(a + z)) and Vl^(aexp(a + 
z)) exp(VF(a exp(a + z))) = aexp(a + ^;) yield 

zir) = r — a + ln(r/a). (7) 

The approximate solution denoted as ipcif) for the spherical system is thus derived as 

, f ^ 1 + 7exp(-(r - a + ln(r/a))) 
i — 7exp(— (r — a + ln(r/ajjj 

where the RHS is the potential in the corresponding plate system with z replaced by z{r). Near a spherical interface 
of large radius, i.e., {r ~ a)/a << 1 and a >> 1, the ipc solution approaches the planar solution, as required. In the 
region far away from the interface (r >> a), Eq. ^ becomes tpci^) = 47acxp(— (r — a))/r, which reduces to the 
linear spherical solution Eq. (|3|) in the weak potential limit. Note that the tpQ solution may be derived algebraically 
by a variable substitution s = a/r exp(— (r — a)) in Eq. ([1]) and more accurate results can be obtained by perturbation 
analysisi^ii^. In comparison to the algebraic method, the derivation of the tpc solution via the geometric mapping 
not only reduces the complexity of algebraic calculations, but also shows how the spherical geometry modifies the 
equipotential surfaces of a planar system as encoded by the geometric deviation d(z). The relation between the ipc 
solution and both the linearized and planar solutions is also revealed in the geometric derivation. 
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Eq. ([8]) is derived from the planar and the hnearized spherical solution, so at the very least it is expected to work 
for either a >> 1 or ^/jq << 1- It is therefore superior to the linear solution, which only works in the weak potential 
limit. The region of validity of the i/'G solution can be derived algebraically. By introducing x = 1/r, the PB equation 
for a spherical system becomes^^ 



where C{x) = x^. Inspired by the functional form of the planar and the linearized spherical solutions Eqs.([3l |4]), we 
use ■(/'(x) = 21n((l + g{x))/{l — g{x))) as a trial solution. Depending on the sign of iJj{x), g{x) = ±exp(— /(a;)) and 
g{x) € (—1, 1) corresponding to V-" G (^oo, oo). Inserting the ansatz into Eq. ^ yields 



For l^l << 1, by dropping the g^ term, the solution to Eq. (|10[) is f{x) ~ l/x — ln.T + ci, with an integration constant 
ci. Inserting f{x) into the ansatz ?/; yields the i/jg solution. An alternative condition for dropping the g'^ term in Eq. 
([T0| is C(a;)(/'^ + /") — 1 << 1, which is equivalent to x = 1/r << 1 by inserting the expressions for f{x) and C{x). 
Therefore, for either \g\ « 1 or x = 1/?' << 1, the solution to Eq. ([T0| coincides with the V'G solution. Note that 
\g\ << 1 is equivalent to the weak potential limit, and x = 1/r << 1 holds for a » 1 since r > a. An important case 
falls in this region of validity of the ipc solution. Consider colloids of size i? in a solution of ion strength / (in mol/L). 
The Debye length is K^^(nm) = 0.304/ /(mol/L) which is at the order of nm for / - 1 mol/lJi. For R - /im, 
a = R/k^^ >> 1. Therefore, the approximate analytical ipa solution is suitable for typical colloidal dispersions. In 
comparison to the series solution^, the formal simplicity of the ipG solution enables further analytical study of the 
electrostatics of colloidal systems. 

The ipG solution turns out to have a larger region of validity. Fig. [3] shows comparisons of the ipc solution (red 
dots), the linearized solution (green dashed) and the numerical solution (black) to the PB equation for a spherical 
system for different potentials. For weak potential {tpo = 0.1(a), 1(6)), the three solutions fall on the same curve. 
The linearized solution works well at least up to 'i/'o = 5 without qualitatively deviating from the numerical solution. 
So the linearized theory applies for moderate values of surface potentiali^. As ipQ exceeds 5, the linearized solution 
starts to deviate from the numerical solution, while Tpc conforms to the numerical solution up to ipQ = 10, where the 
linearized solution deviates significantly the numerical solution. This indicates that the dependence of d{z) on -00 is 
weak for a = 1 up to tJjq = 10. 

The quality of the "ipG solution can be systematically studied by defining a ratio 6 = max^lKA-!/' — smhip)/ smhip\}. 



C{x)—^iP = smhij, 



(9) 



ff(x)[C(x)(/'2 - /") - 1] + g'[C{x)ir + ./") - 1] = 0. 



(10) 




FIG. 4: The ipc solution is applicable below the red curve in the parameter space {a, -^o}- On the red curve, 5 = e = 10 ^. 

The smaller the ratio S is, the better the solution is. For a given precision e = 10^'^, the applicable region of the 
t/'g solution is found to be below the red curve in the parameter space {a, V'o}: as shown in Fig. |4l For a >7, the 
ipG solution applies even for large potentials. There exists, however, a cut-off value for the surface potential. High 
potential, or equivalently low temperature, may lead to correlation of counter-ions near the charged interface that is 
ignored in the mean field PB equation^. In addition, high potential leads to high concentration of ions so that the 
finite dimension of ions must be taken into consideration!^. The advantage of the tpc solution over the planar solution 
is shown explicitly in Fig. [Sj The ipc solution works better than the planar solution even for large spherical radius. 
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FIG. 5: The comparison of the tpo solution (black) and the planar solution (blue) in terms of the ratio of the LHS to RHS of 
Eq. The radii of the spherical interface are a — 5(a) and 20(6). tpo = 1. 

In conclusion, we have studied the EDL structure around charged spherical interfaces by analysis of the Poisson- 
Boltzmann equation. Despite the point charge assumption of electrolyte ions and the neglect of ion-ion corrections, 
the PB equation generally works well especially for problems of electrostatic interaction of colloidal particles^. In 
this paper, we have derived an approximate analytical solution to the Poisson-Boltzmann equation for the spherical 
system by a geometric mapping. The formal simplicity of the ^pa solution enables further analytical study of spherical 
systems. The regime of applicability includes not only the weak potential regime where the linearized solution also 
works well, but also the regime of large spherical radius. Typical colloidal dispersions with the size of colloids much 
bigger than the Debye length fall in the latter regime. 
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